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SYMBOLS 


blowing parameter defined by equation [19) 
heat transfer coefficient 
specific heat 
total enthalpy 

convection coefficient parameter 

volumetric convective heat transfer coefficient 

external incident radiation intensity 

radiation intensity outward in the matrix 

radiation intensity transmitted inward in the matrix 

absorption coefficient 

thermal conductivity 

thickness of the matrix 

molecular weight 

index on radiation bands that penetrate the matrix 
mass injection rate 

index on radiation bands absorbed on the exposed surface 
pressure 

normalized pressure 
heat transfer rate 
universal gas constant 
scattering coefficient 
temperature 
time 

velocity in x direction 
distance from front surface 

iii 
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viscous coefficient 
inertial coefficient 
emissivity 
relative error 
reflectivity 

cooling effectiveness defined in equation (42) 

normalized temperature 

viscosity 

normalized spatial coordinate 
density 

Stefan-Boltzmann constant 
porosity 

convective heat transfer correction factor by mass injection (see 
eq. (17)) 

radiative heat transfer correction factor by mass injection (see 
eq. (20)) 


Subscripts : 
b back face 

aond conductive 

e external 

e edge of boundary layer 

f fluid 

i internal 

o without injection 

p radiative 

radin incident radiative 

verad reradiative 


IV 


s 


solid 


w wall 

c» free stream 

V the vth spectral band of radiation 

Superscripts : 

( ) ' differentiation with respect to 5 
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A SIMPLIFIED ANALYTICAL SOLUTION FOR THERMAL RESPONSE OF A 


ONE-DIMENSIONAL, STEADY-STATE TRANSPIRATION COOLING 
SYSTEM IN RADIATIVE AND CONVECTIVE ENVIRONMENT 
Hirotoshi Kubota* 

Ames Research Center 


SUMMARY 


A simplified analytical method for calculation of thermal response within 
a transpiration-cooled porous heat shield material in an intense radiative- 
convective heating environment is presented. The essential assumptions of the 
radiative and convective transfer processes in the heat shield matrix are the 
"two-temperature" approximation and the specified radiative-convective heat- 
ings of the front surface. Sample calculations for porous silica with CO 2 
injection are presented for some typical parameters of mass injection rate, 
porosity, and material thickness. The effect of these parameters on the 
cooling system is discussed. 


INTRODUCTION 


Probes entering the atmospheres of large planets, such as Jupiter and 
Saturn, encounter severe radiative heat loads in addition to the usual convec- 
tive heat loads. The ablating heat shields, both charring and subliming types, 
used in the past to protect probes from convective heating, may not be ade- 
quate for these additional heat loads. For such entries, ablators will have 
to be chosen to withstand both the radiative and the convective heating. 

One method, which is now being developed to withstand the radiative and 
convective heat load, uses reflecting ablation heat shields [refs. 1-3). It 
should be noted, however, that certain ablators, for example, silica, which 
are effective in protecting against convective heat may not block radiative 
heat and may, in fact, increase the incident radiative heating at certain 
ablation rates (ref. 4} . 

Another method which can be easily adapted to withstand the radiative and 
convective heat loads is a transpiration cooling technique. By choosing a 
porous heat shield which has a high reflectivity, part of the radiative heat 
load can be reflected. By choosing an absorbing transpiration gas (ref. 5), 
part of the radiative heat load can be absorbed and the convective heat load 
can be blocked. 


*National Research Council Associate at NASA Ames Research Center, 
Moffett Field, California, 94035, during the author's official leave from 
National Aerospace Laboratory of Japan. 
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Accurate predictions of heat and mass transfer in the porous media must 
be solved from the coupling of momentum and energy equations for fluid and 
energy equation for a solid matrix (ref. 6). Since the convective heat 
exchange between solid and fluid affects the temperature distributions of solid 
and fluid in the porous matrix, two basic techniques have been introduced. One 
is the "one-temperature" assumption, that solid and fluid temperatures are the 
same, for a large heat transfer or small mass transfer rate, such as is the 
case in charring ablation (refs. 7 and 8); the other is the "two-temperature" 
assumption, that solid and fluid temperatures are not the same, an assumption 
that is essential for a small heat transfer and a large mass transfer rate 
(ref. 9) . 

In this investigation, the transpiration cooling system is studied 
theoretically with the intent of finding simple analytical solutions. This 
involves the interactions of the heat and mass transfer process within the 
porous heat shield. For this purpose, a simplified, one-dimensional, steady- 
state case with constant properties will be investigated. It is expected that, 
from this simplified model, the essential features and the basic parameters 
that govern the transpiration cooling system can be found. 


ANALYSIS 


Physical Model and Assumptions 

The physical model considered in this analysis is sketched in figure 1. 

A one-dimensional porous matrix of thickness L is exposed to the external con- 
vective and radiative heat fluxes. The present analysis considers the heat 
transfer process only within the porous matrix under the conditions of given 
convective and radiative heating at the front surface. The coolant gas is 
injected through the matrix from back to front surfaces. 

The analysis makes the following major assumptions: 

1. Steady-state transpiration 

2. Solid and gas temperatures in the matrix not equal (two-temperature 
assumption] 

3. Continuum model for the porous matrix (volume element is large with 
respect to individual pores) 

4. No reversal flow (all injected gases flow toward the heated front 
surface) 

5. Specified incident heating fluxes and constant back-face temperature 

6. Constant thermal and optical properties (independent of temperature, 
pressure, and wavelength) 
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7. Neither chemical reactions between the porous matrix and the injected 
gas nor dissociation of the injected gas 

8. Negligibly small energy change due to the pressure difference in the 
porous matrix in comparison with other energy terms 

9. Front surface impermeable to the external gas 


Simplified Governing Equations 

The complete transient differential equations, describing the conserva- 
tions of mass, momentum, and energy without radiation in the porous matrix and 
with mass transfer, are presented in reference 10. By adding a radiative 
transfer term to the energy equation, the differential equations can be 
applied to the case with radiation. This radiative transfer term will be 
derived in the next section. The complete equations will not be presented 
here. The simplified equations, obtained by applying the previous assumptions, 
are as follows. 

For fluid: 

Mass conservation 


dm 

dx 


= 0 


Momentum conservation - Darcy's equation 




Energy conservation 


d^T 


dT, 


*=/ " V I - V - nd 




Perfect gas law 






For solid: 
Energy conservation 


d'^1 a fi 






( 1 ) 


( 2 ) 


(3) 


(4) 


( 5 ) 


where a, B, and j) are constants that depend on the porosity, and ih is a 
volumetric heat transfer coefficient between solid and fluid and can be 
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written by the experimental correlation (ref. 11) as 




1 .44 


( 6 ) 


From equation (1) one obtains 

m = ~Pf^f ~ constant 

that represents a given coolant injection rate at the back face. From equa- 
tions (2) , (4) , and (7) , one also obtains 


(7) 


dp 


'f 


Py. = id\if + ^m')m{R/Mj^)T 

Then, equations (3), (5), and (8) become a set of the governing equations. 


( 8 ) 


Radiat-ive transfer term, dq^/dx.- The matrix material is partially trans- 
parent to radiation so that a part of the incident radiative flux penetrates 
the material to be backscattered in depth and the rest is absorbed at the sur- 
face (see fig. 2). The incident radiation is defined (ref. 2) as 


m' 


n ' 




(9) 


v=m 


v=n 


where the first sum is radiation that penetrates and the last is surface 
absorbed radiation. 

Thus, the radiative flux within the matrix should ^e evaluated under 

the influence of the penetrating incident radiative flux qy,^ ; hence, by 

v=m ^ 

the "two-flux" method described in appendix A, it is written as 


m' 


m' 


= T. ^r.= T, ■ V 


( 10 ) 


v=m 


v=m 


The radiative terms in equations (3) and (5) are then expressed, from the 
solution of the Kubelka-Munk equations (Al) and (A2) , by 


dPy, 

dx 


m' /dim dip \ m' 

v=m ^ ' v-m 


+ dp ) 

V "v 


( 11 ) 


(^) 


m' 

-2 \ ^]e^K ^r^ ^ ^ 

v=m V 




m’ 


g-axp qy. 


v-m 


( 12 ) 


4 



Boundary conditions.- Boundary conditions required for the simplified 


governing equations 

are (see fig 

■ 3) 


at X = 0; 

dTc. 

dx 

~ ^radin ~ ^rerad 

(13) 


d^T.p 

kf ^ 

^ dx2 

dT. 

- 0 or -kj^ dx ~ ^ 

(14) 


Pf 

II 

(15) 

at X = L', 

Ts 

II 

11 

(16) 


where the convective heating rate with blowing (ref. 12), the incident radia- 
tive heating flux with' blowing^ and the surface reradiation heat flux are 
defined as 


= Pe^e^Ho'^con^^e ' 


'^WQ^^aon 




con B 


= 


e ^ - 1 


m 




dradin £ (1 " ^rad 

v=n V 


rerad 


- eaTc 


■'w 


(17) 

(18) 

(19) 

( 20 ) 
( 21 ) 


In equations (17) and (20), ^^^n '^rad correction factors for 

incident convective and radiative fluxes by mass injection, respectively. 


Nondimensional Governing Equations and Boundary Conditions 

The simplified equations (3) , (5) , and (8) can be written in the 
nondimensional forms as follows: 

+ 3>5 = 0 ( 22 ) 

^Private communication with J. T. Howe. 
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where 


0^ + $20} + ^3C0s - 0/) + ^6 + ^7 = 0 

p^Py. - $ 80 / - ®9 = 0 


. _ Ts-^b 

, _ ^ f - '^b 

^f ~ T - Tu 


L 


Pf Pf 


w 




(23) 

(24) 

(25) 

(26) 

(27) 

(28) 


$1 = -hL^lkgH - <(.) 
«>2 = 

$3 = hL^lk^ 


m' 


$c = 


2KL'^ {Aq/Ai) y] q-p^ 
\>=m ^7 

- ( 4 iSo/ 4 oBi )$4 




$6 = (?Cg/?Cy )$4 


0 7 


On — 


- {kg/k^) {A\Bq/ AqBi)^i^ 
Lm{B/Mf){Tg^ - 7’^)(dy^ + Bm)}p^ 


(29) 


( ) ' : differentiation with respect to C 
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The corresponding boundary conditions (eqs. (13) through (16)) are 


at C = 0; 

6s (0) = 1 

(30) 

where 

08 ^Tadin ~ ’^revad) 

(31) 


$10 = -L/k^iTg^ - %-) 

(32) 


CD 

0 

II 

0 

(33) 

or 

CD 

0 

II 

0 

(34) 


1— 1 
II 

0 

(35) 

at C = 1] 

0^(1) = 0 

(36) 


0 

II 

r-H 

CD 

(37) 

In equation (31), 
temperature, Tq^. 

Evadin’ ^rerad functions of only the 

Method of Solution 

surface 


General solutions of equations (22) and (23) can be written as 

+ E2 


(i: 

\k=i 


08 " I ^ 


) 


(38) 




Kk=i 


) 


0f(?) = ( H -Tk^k e e 


-aL^ 


(39) 


where the bracketed terms are the homogeneous solutions and the other terms 
are the particular solutions. 


Coefficient 6^ in equations (38) and (39) is obtained from the homogene- 
ous solutions of equations (22) and (23) as discussed in appendix B. The 
other coefficients are obtained from the boundary conditions, equations (30) 
through (37). Note that all coefficients except 6^ are implicit functions of 
7'g^^. The value is determined by iteration to satisfy the energy balance 

given by equation (31) and equation (38) for ^ = 0, that is. 


3 

0g'(O) = ^2 ^k^Jc * “ £’ 2 “-^ (40) 

k=l 
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The pressure variable is solved from equation (24) by integrating 


p/(o 


“ 3 

2$8 Z ^ Cff3/aL)(e“^^ - l) 

- (^^/aL) - l) + 2{y^c^^Q + $9)? + 1 


Sample Calculations 

The effect of radiation penetration into the porous matrix has been 
tested for two boundary conditions, corresponding to (1) no conduction heat of 
the fluid at the front surface and (2) zero derivative of the fluid conduction 
heat at the front surface. The following table shows four typical sample 
calculations considered in this analysis. 


Case 

Radiative penetration 

Boundary condition 

I-A 

Without 

e^”(0) = 0 

I-B 

Without 

Bf' (0) = 0 

II-A 

With 

CD 

0 

II 

0 

II-B 

With 

0^' (0) = 0 


r 


In this sample calculation, silica (Si02) and carbon dioxide (CO 2 ) are 
taken as porous matrix and transpirant gas, respectively. Detailed descrip- 
tion of their thermal and optical properties are shown in appendix C. 

The radiative and convective heat loads for this sample calculation are 
taken for a 15° Saturn entry (ref. 3) at t = 23.7 sec, corresponding to maxi- 
mum radiative heat loads. 

= 216.70 cal/sec/cm^ 

^ o 

p = 2.091 atm 
^ e 

U = 2.86x10® cm/sec 

00 

PgU^Cf] = 1.186x10“^ gm/cm^/sec 
Hg = 0.9685x10® cal/gm 
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The spectral radiative heating data for Si02 blowing can be broken into 
20 groups whose wavelength ranges cover 0.0919 a through 9.92 p. (The exact 
incident heating for CO 2 injection may be also estimated by the use of some 
numerical computations such as RASLE CODE (ref. 4).) The incident radiations 
are combined into groups 1-3 (a part of infrared) and 17-20 (ultraviolet) that 
are absorbed at the surface of the silica; groups 4-16 penetrate into the 
matrix. ^ 


RESULTS AND DISCUSSIONS 

Typical sample solutions are presented for a mass injection rate of 
w = 0.1 gm/cm^/sec, ~ = 0- The solid and fluid tempera- 

ture distributions and the pressure distributions for these sample cases are 
shown in figure 4(a) through figure 4(d) for the porosity range from (j) = 0.2 
to (j) = 0.8. As would be expected, temperature distributions Tg and Ty? are 
sufficiently close to each other so that the "one- temperature" approximation 
will provide a good solution at the smaller porosity. However, the "two- 
temperature" method is required for higher porosity cases. 

The effect of radiative transfer within the porous matrix on temperature 
and fluid pressure is presented in figure 5(a) for the same input and boundary 
condition as before, except for a porosity of ()> = 0.4. A significant differ- 
ence in temperature between two cases, with and without radiation heat, is 
observed . 

The effect of boundary condition (case A, 0y"(O) = 0; and case B, 

(0) = 0) on solid and fluid temperature distribution is shown in figure 5(b) 

for the same input conditions as in figure 5(a) except for a larger mass injec- 
tion rate of m = 0.2 gm/cm^/sec. There is no significant effect of the bound- 
ary condition on the solution except in the vicinity of the front surface. 

Since the boundary condition, 6f"(0) = 0, is a more reasonable assumption 
(ref. 10), the remainder of the calculations will be made under this 
condition. 

Figure 6 presents the cross-plot of wall temperatures and back-face 
pressure of fluid as a function of porosity for various radiation heat trans- 
fer correction factors >p 2 ^ad- Variations of the blowing correction factor for 
the radiative heating 'CpatZ with blowing parameter is given for Si 02 and carbon 
ablation (ref. 4) . 

The significant reduction of fluid wall temperature due to the matrix 
porosity is noted, as well as larger separation of solid and fluid tempera- 
tures as the value of porosity increases. 

The variations of solid surface temperature due to mass injection 

rate and porosity, are shown in figure 7(a). Note that Tg is not signifi- 

h) 

cantly affected by the change of porosity to a higher value (ej) ^ 0.4). This 
effect can be also recognized partly from figure 6. 

^Private communication with J. T. Howe. 
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The cooling effectiveness is demonstrated in figure 7(b) for the various 
porosity and mass injection rate. The cooling effectiveness is defined by 


n 



(42) 


where Tj^, and Bq are wall temperature with and without mass injection, 

back-face temperature, and injection parameter Bq = m/ respectively. 

The cooling effectiveness n depends strongly on the mass injection rate, 
decreasing as the injection rate increases. However, the effect of porosity 
on cooling effectiveness, for a given mass injection rate, is not so signifi- 
cant; that is, the cooling effectiveness may be expressed in a single average 
curve (e.g., (]) = 0.4 case) that represents overall correlation between solid 
temperature, porosity, and mass injection rate. 

The effect of material thickness on wall temperature and back-face heat 
conduction is shown in figure 8. The material thickness has an important 
effect on reducing back-face heat conduction but it has practically no effect 
on solid surface temperature. 

The characteristics of the solutions, described above, can be plainly 
explained by using the concept of the convection coefficient parameter 
(ref. 13), H, defined by 


H = 



(43) 


where h, the heat transfer coefficient, is inversely proportional to the 
matrix porosity. The parameter H, .given by equation (43) , signifies the 
degree of the heat exchange between the solid and fluid; that is, the larger 
the convection coefficient parameter (larger heat transfer coefficient or 
smaller porosity; larger matrix thickness or smaller mass injection rate, or 
both) the less difference there is between the solid and fluid temperatures. 


The important criteria for a transpiration cooling system are cross- 
plotted in figure 9 with the assumption that the porous silica used for the 
present cooling system will melt at 1981 K (ref. 3). Then, the boundary 
criteria, in terms of radiation heat transfer correction factor 

porosity limiting temperature, can be cross-plotted as in figure 9. In this 
figure, the region above each boundary curve requires ablation cooling while 
the other region requires only transpiration cooling. Also shown in the same 
figure is the functional curve, the vs m relation, obtainable from the 

external flow condition (ref. 4) . 

The cross points between the solid lines and the dashed line represent 
the minimum mass injection rate required for maintaining a porous material 
below the melting temperature. The higher porosity requires less mass injec- 
tion rate for the given material thickness. 
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The thermal properties have important effects on the solutions. 

Figure 10 is a composite plot showing the effects of porosity, material thick- 
ness, and thermal properties on the front-surface solid and fluid temperatures. 
Nondimensionalized temperature differences for the solid and the fluid at the 
front surface, due to thermal properties (see appendix C) , are defined as 


Tg^Ccase <a>) - Tg^Ccase <b>) 




^SyCcase <a» 




w 


|rj?^(case <a» - Tjy^(ca.se <b>) 

Tf (case <a^>) 

•‘w 


(44) 


respectively. 


Note that the thermal properties for case <a> are selected at the melting 
temperature of silica, that is, 1981 K. For case <b>, the back-face tempera- 
ture of 300 K is chosen. These are two extreme temperatures expected. Since 
zrp and are the relative errors in temperature as a result of the 

T 1 } 

thermal properties, the magnitude of these parameters shows the effect of the 
thermal properties. Note again, that the maximum effect is presented because 
the two extreme temperatures were taken. 


The effect of the thermal properties is quite important on the solid wall 
temperature but not too important for the fluid wall temperature. The effect 
of porosity is to increase the temperature difference. The thinner material 
thickness is also strongly affected. 


Finally, the present analytical solution has been compared with a more 
rigorous numerical solution (but still with constant thermal and optical prop- 
erties) that retains the energy change term due to the pressure difference in 
the porous matrix (see assumption 8) . The result shows that the present 
analysis compares within 1 percent of the rigorous solution, that is, the 
Adams-Moulton integral scheme. This result verifies the assumption that the 
term Uj^[dpj^/dx) can be neglected except for the case of extremely small 
porosity, that is, larger pressure drop within the porous matrix. 


CONCLUDING REMARKS 


A simplified analytical method for calculation of thermal response within 
a transpiration-cooled heat shield material in an intense radiative-convective 
heating environment is presented. 

Essential features of the analysis are the "two- temperature" assumption, 
where matrix and flowing gas have different temperature distributions, and the 
specified radiative-convective heatings at the front surface. 
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Sample calculations for porous silica with CO 2 injection result in the 
following major conclusions: 

1. Radiative transfer within the heat shield material increases the 

temperature of the material. Therefore, the absorption of radiation by trans- 
pirant gas does not play an effective role as a heat shield, but, the injec- 
tion of the gas into the adjacent boundary layer reduces the incident heating 
and the temperature level within the matrix becomes lower (effect of and 

'^rad) • 

2 . Increasing the mass injection rate of the transpirant gas reduces 
temperature distributions, but increases the fluid pressure difference between 
front surface and back face. Cooling effectiveness per mass injection rate 
decreases as mass injection rate increases. In this respect, the behavior is 
similar to that of the general mass transfer cooling system. 

3. Small porosity increases the volumetric heat transfer coefficient, 
thus causing closer temperature distributions between solid and fluid and 
resulting in higher back-face fluid pressure. 

4. Increasing material thickness gives lower back-face heat conduction 
but has practically no effect on solid surface temperature. 

5. The thermal conductivity and specific heat of the solid matrix have 
more important effects on the solutions than fluid thermal properties. 

Agreement of the present solution with the more rigorous numerical 
calculation (the Adams-Moulton integral scheme) is quite good. 

The calculated results were presented in the terms of the correlational 
form for mass injection rate, porosity, and radiation heat transfer correction 
factor, as the criteria for a transpiration cooling system. 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, September 16, 1975 
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APPENDIX A 


RADIATIVE TRANSFER EQUATIONS 


When the radiative emission inside the matrix is negligibly small 
compared to the radiant fluxes transmitted inward and backscattered outward, 
the one-dimensional radiation transfer can be described by the Kubelka-Munk 
"two-flux" equations (ref. 14) as 

djy 

(Al) 

diR 

^ = (5 . - SlT^ (A2) 


where K and S are absorption coefficient and scattering coefficient, 
respectively, and are defined by 


K = (1 - <t>)Kg (A3) 

S = ^ (1 - <|,)5s (A4) 

The boundary conditions for this set of simultaneous equations (Al) and 
(A2) are 

at r = 0 (front surface) ; 

at X = L (back face) ; 

as illustrated in figure 2. 

The general solutions are 


= 4(1 - 6)e“^ + S(1 + 3)e (A7) 

Ir^ = 4(1 + 3)e“^ + B(1 - 3)e"“^ (A8) 
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I II II 


l■ll■■lll■ll mill 




mill II 


where 

a = ^K{K + 2 S) 

3 = ^/K/{K + 2 S) 

^ 0 

A = (1 - 

Bo 

B = -(1 - 

-olL 

^0 = [1 - s - 02,(1 + S)]e 

Ai = [(1 - B)2 - (1 - + 02,) + (1 + B)^02:02,]e 

. 9 1 QlLj 

-[Cl + 3)^ - (1 - 3^) (o^ + 02,) + (1 - 3) o^02,]e 

Sq = (1 - 3^) Cl + 02,^) " 2C1 + 3^)02, 

B\ = ^ 0^1 
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APPENDIX B 


EXPRESSIONS FOR COEFFICIENTS IN SOLUTION 


Characteristic Equation and Its Solution 
By substituting the bracketed terms of equations (38) and (39) into the 


homogeneous forms of equations (22) and (23) , one obtains 



+ $l(c^ - = 0 

ik = 2, . . 

k') 

(Bl) 

-^k^k^k * ^^-^k^k^k ^ ^^^^k - ^k^k^ = 0 

{k = 1 , 2, . . 

k') 

(B2) 

Yo = 1 



(B3) 

Then, since Ci^ f Q , 




Yfe = ^ 4>l)/«>l 

{k = 1,2, . 

. k') 

(B4) 


6^3 + $25^^ + ($1 - $3)<S?. + = 0 (k= 1,2,. , k') (B5) 


and k' =3. Therefore, 6^ are the roots of the characteristic equation 

63 + $262 + = 0 (B6) 


Coefficients of Solutions for Case 1-A and Case II-A 

By substituting equations (38) and (39) into equations (22) and (23) , one 
obtains 


El = -(a2$6 + + CCiai^) 

E2 ~ - (i>2^7 ^4’^5) / ('^2'^ 3 ^1^4) 

£"3 = (a3«>4 - ai$6)/(a2a3 + aian) 
S’4 = (&3$5 - &l^’7)/(^2^3 ^1^4) 


(B7) 


where 
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a\ = + $1 

a-i ~ 

as = $3 

(24 = + ^'loL - $3 

hi = + $1 


&2 = 
bs = 3>3 

2)4 = - $2'^'^ - 5*3 j 


(B8) 


By applying the boundary conditions, equations (30), (33), (36), and (37), 
to equations (38) and (39) 


where 


^■1 

A2 

^43 

/ 


F2 

Bl 

Bl 

BET 

(B9) 

Ps 

Cl 

Cl 

/ 



Fl 

^1 

/ 


Bl 

Fl 

Bl 

BET 

(BIO) 

Cl 

Fl 

Cl 

/ 


'.Ai 

^2 

Fl 

/ 


Bl 

Bl 

Fl 

BET 

(BID 

Cl 

Cl 

Fl 

/ 



^0 


1 - El 

- E2 - 

■ 

- C2 - Os 
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For case II-A: 


For case I-A: 


= 1 - e ^ ik = I, 2, 5) (B14) 

St, 

= 1 - e (fe = 1, 2, 3) (B15) 

Ck = (fe = 1, 2, 3) (B16) 

F’l = 1 - fflCl - e“^) - £"2 Cl - CB17) 

F2 = 1 - El - E2 + Es e“^ + £"4 CB18) 

£-3 =- a 2 £ 2(^3 + 

£■]: = £"2 = ffs = £4 = 0 (B20) 


Coefficients of Solutions for Case I-B and Case II-B 

The boundary condition (33) of case I-A and case II-A is replaced by 
equation (34) , therefore equations (B16) and (B19) should be replaced by 


= yj^^T^ {k = 1, 2, 3) (B21) 

£3 = aL{-E^ + £ 4 ) (B22) 


respectively. 
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APPENDIX C 


THERMAL AND OPTICAL PROPERTIES OF POROUS MATRIX 
AND TRANSPIRANT GAS 


Thermal and optical properties of porous silica and carbon dioxide used 
the matrix and transpirant gas are given in tables 1 and 2. 


TABLE 1 - PROPERTIES OF POROUS SILICA 





Value at 


Reference 

Property 

1981 K 
(case <a» 

300 K 
(case <b>) 

Dimension 

kg 

5.470x10"3 

2.157x10”3 

cal/sec/ cm/K 

3 

^Ps 

0.309 

0.172 

cal/gm/K 

3 

P 

1.49 

1.49 

gm/cm^ 

3 

s 

0.001 

0.001 

cm" ^ 

3 

Ss 

40 

40 

cm"^ 

3 


0.06511 

0.06511 


3 

hi 

0.17 

0.17 


3 

'^b 

0.8 

0.8 


3 

e 

0.14 

0.14 
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TABLE 2 - PROPERTIES OF CARBON DIOXIDE 


Property 



Value at 

Dimension 

Reference 

1981 K 
(case <a» 

300 K 
(case <b» 


3.482x10”'* 

2.387x10-5 

cal/ sec/cm/K 

12 


0.327 

0.202 

cal/ gm/K 

16 

Kf 

0.1172 

0.1172 

cm” ^ 

17 


0 

0 

cm” ^ 

17 


5.989x10”'+ 

1.500x10”*+ 

gm/ cm/ sec 

12 





Viscous Coefficient and Inertial Coefficient 


Since the data for viscous coefficient 5 and inertial coefficient 3 
not available for porous silica, the experimental correlations for porous 
feltmetal (ref. 9) 


a = 0.7336x10^ (100 (}>) 1- 5378 (cm’^) 
3 = 0.1517x10^+ (100 (j))0-2520 (cm"l) 


which are obtained from the extrapolated data, are used. 
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Figure 2.- Illustration of radiative transfer model. 
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Figure 5.- Surface energy balance. 


24 



Front surface 


Back face 
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Figure 4.- Temperature and fluid pressure distributions within porous matrix: 
m = 0.1 gm/cm^/sec, L = 1.0 cm, ip , = 1.0, case II-A. 
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(a) Effect of radiative penetration for case A: 

m = O.l gm/cm^/sec, L = 1.0 cm, ^ = 0.4, 

ij; , = 1.0. 

rad 


(b) Effect of boundary condition (with radiative 
penetration): m = 0.2 gm/cm^/sec, 

L = 1.0 cm, (f) = 0.4, - l-O. 


Figure 5.- Temperature and pressure distributions. 





Figure 6.- Solid and fluid temperature variation at wall and fluid pressure at 
back face as a function of porosity: m = 0.1 gm/cm^/sec, L = 1.0 cm, 

case II-A. 
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Figure 8.- Effect of material thickness; solid 
surface temperature and back-face heat con- 
duction: (j) = 0.4, ^ ~ l-O? case II-A. 
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Figure 9.- Criteria for a transpira 
system (desired minimum mass inj 
case II-A, L = 1.0 cm. 
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